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We derive an upper limit of Bo < 3.4 x 10“®(Oo^5o)^ Gauss on the present strength of any primor¬ 
dial homogeneous magnetic field. The microwave background anisotropy created by cosmological 
magnetic fields is calculated in the most general flat and open anisotropic cosmologies contain¬ 
ing expansion-rate and 3-curvature anisotropies. Our limit is derived from a statistical analysis of 
the 4-year Cosmic Background Explorer data for anisotropy patterns characteristic of homogeneous 
anisotropy averaged over all possible sky orientations with respect to the COBE receiver. The limits 
we obtain are considerably stronger than those imposed by primordial nucleosynthesis and ensure 
that other magnetic field effects on the microwave background structure are unobservably small. 
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The origin of magnetic fields observed in galaxies and 
galaxy clusters is still a mystery. The invocation of pro- 
togalactic dynamos to explain the magnitude of the field 
involves many uncertain assumptions but still requires a 
small primordial (pregalactic) seed field ||]. Hence the 
possibility of a primordial field merits serious consider¬ 
ation. Other attempts to find an origin for the field in 
the early universe have appealed to battery effects, the 
electroweak phase transition, or to fundamental changes 
in the nature of the electromagnetic interaction. All in¬ 
troduce further hypotheses about the early universe or 
the structure of the electroweak interaction [|j. All aim 
to generate fields by causal processes when the Universe 
is of finite age. Therefore, any magnetic field created by 
these means will exist only on very small scales with an 
energy density that is a negligible fraction of the back¬ 
ground equilibrium radiation energy density. 

Nevertheless, while such fields might still provide the 
seeds for non-linear dynamos in the post-recombination 
era, any large-scale primordial magnetic field with a 
strength of order B ~ 10“® Gauss, comparable to that 
inferred from the lowest measured intergalactic fields and 
close to the observational upper limits via Faraday rota¬ 
tion measurements |^, may well be of cosmological ori¬ 
gin. A similar pregalactic (or protogalactic) field strength 
is inferred from the detection of fields of order 10“® Gauss 
in high redshift galaxies and in damped Lyman al¬ 
pha clouds Q, where the observed fields are likely to 
have been adiabatically amplified during protogalactic 
collapse. In the absence of a plausible dynamo for gen¬ 
erating large-scale pregalactic fields, it is of interest to 
reconsider the limits on a large-scale primordial field in 
view of new observational constraints that we outline be¬ 
low. 

Primordial magnetic fields can leave observable traces 
of their influence on the expansion dynamics of the Uni¬ 
verse because they create anisotropic pressures and these 
pressures require an anisotropic gravitational field to sup¬ 
port them. Primordial nucleosynthesis constraints only 
limit the equivalent current epoch field to be less than 
about 3 X 10“^ Gauss 0, a value that is only slightly 


stronger than the dynamical constraint at nucleosynthe¬ 
sis ^ . We show in this letter that the cosmic microwave 
background isotropy provides a stronger limit on the 
strength of a homogeneous component of a primordial 
magnetic field. 

We consider the cosmological evolution of the most 
general homogeneous magnetic fields, calculate their 
gravitational effects on the temperature anisotropy of 
the microwave background radiation, and hence derive 
a strong limit on the strength of any homogeneous cos¬ 
mological magnetic field by using the 4-year Cosmic 
Background Explorer (COBE) microwave background 
isotropy measurements ||^. We employ statistical sam¬ 
pling techniques appropriate for the non-Gaussian statis¬ 
tics of the large-scale temperature anisotropy pattern cre¬ 
ated by a general homogeneous cosmological magnetic 
field and allow for the randomness of the angle at which 
COBE views the characteristic anisotropy pattern on the 
sky 0 - The addition of a homogeneous cosmological 
electric field will not be considered: a homogeneous in¬ 
tergalactic electric field would create a current of charged 
particles, and rapidly decay. 

Our limit derives from the non-linear coupled evolu¬ 
tion of the shear anisotropy and magnetic field density 
during the radiation era which we shall discuss below. 
In the presence of an equilibrium background of black- 
body radiation with isotropic momentum distribution, 
pressure Pr, density pr, and equation of state Pr = ^Pr, 
the anisotropic magnetic pressure prevents the rapid de¬ 
cay of the expansion shear anisotropy familiar in Bianchi 
type I universes containing only matter with isotropic 
pressure. An approximate analysis of the problem in an 
axisymmetric magnetic Bianchi type I universe was given 
by Zeldovich j^. The isotropic expansion is stable at 
second-order and the anisotropies decay only logarithmi¬ 
cally in time relative to the mean expansion rate. Phys¬ 
ically, the evolution of the anisotropy is governed by the 
magnetic pressure anisotropy. 

We need to describe the evolution of the magnetic field 
strength and the accompanying shear anisotropy, which 
will distort the microwave background isotropy, in the 


1 



most general anisotropic universes possessing shear and 
3-curvature anisotropies. From earlier studies of the evo¬ 
lution of the most general homogeneous universes in the 
absence of magnetic fields the pattern of evo¬ 

lution of magnetic universes during the radiation era 
can be deduced. In general, if the universe contains a 
non-interacting mixture of isotropic blackbody radiation 
and any matter source possessing an energy-momentum 
tensor with zero trace (e.g. magnetic or electric fields, 
long-wavelength gravitational waves, or an anisotropic 
distribution of collisionless massless or relativistic par¬ 
ticles ll^), then the evolution follows the same char¬ 
acteristic pattern. The most general anisotropic flat 
and open universes containing isotropic universes are of 
Bianchi type VII and are equivalent to simpler flat or 
open anisotropic spacetimes of Bianchi types I or V to 
which gravitational waves have been added. The Ein¬ 
stein tensor for the more general models can be split 
into two pieces: one corresponding to a simpler Bianchi 
type universe, the other to an effective energy-momentum 
tensor for the gravitational waves which has zero trace 
|^,Q. This means that the anisotropies in simple 
Bianchi type I (flat) or type V (open) universes contain¬ 
ing blackbody radiation plus trace-free matter sources 
with anisotropic pressures have similar time-evolution to 
the most general anisotropic universes of type VII con¬ 
taining blackbody radiation and a magnetic field: the 3- 
curvature anisotropies and the magnetic stresses display 
the same asymptotic time evolution. When the devia¬ 
tions from isotropy are small, the evolution of the shear 
and the energy density of any trace-free matter field with 
anisotropic pressures is well-approximated by setting the 
blackbody density and the volume-averaged Hubble ex¬ 
pansion rate equal to their values in the isotropic flat 
Friedmann universe {pr = and H = since we set 
SttG = 1) before solving for the non-linear evolution of 
the shear anisotropy and the magnetic field. The ra¬ 
tio of the shear to Hubble expansion rate has a generic 
behaviour for p < ip; (when p > the anisotropic 
stresses dominate at large t and the solution ceases to be 
a small perturbation of an isotropic Friedmann universe 
as t —*■ oo). In general, when anisotropies are small, the 
ratio of the shear anisotropy, tr, to the mean Hubble rate, 
H, relaxes to a constant value determined by the ratio of 
the total energy densities in anisotropic traceless fluids 
and 3-curvature anisotropies, Pga, and magnetic fields, 
Pb, to that of the isotropic perfect fluid density, p, (radi¬ 
ation or dust). When the isotropic fluid has equation of 
state p = (7 — l)p, with 0 < 7 < 4/3, the time evolution 
is determined by the Einstein equations 



Here, ps = is the magnetic field density. Hence, 


as t ^ 00 , we have a/H 

/ Pga + Pb \ _ / 

dt \ p ) \ 



—> constant and 
'' Pga + PB ^ 

16 ^ ^ Pga+Ps '^ 
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In the radiation era (p = pr,^ = 4/3), at redshifts 1 + z > 
1 + Zeq = ‘&y~ I0^Ho/i|q, where Hq < I is the cosmological 
density parameter and /150 is the present value of the 
Hubble constant in units of 50Kms~^Mpc~^, we have 


Ps/Pr ^ Q/[l + 4(5ln(t/to)]; Q constant. (3) 


During the dust era (p = pd ,7 = 1), when z < Zeq, 
the evolution is determined at linear order and {pga + 
Pb)/P d oc 1 - 1-2 falls linearly with redshift. In general, 
the shear distortion created by magnetic fields and any 
other trace-free anisotropic stresses is given by 
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where <5 is a constant. The 5 term gives the simple shear 
decay for universes with isotropic 3-curvature contain¬ 
ing only matter with isotropic pressure; this term be¬ 
comes negligible at late times. In both dust and ra¬ 
diation eras the anisotropic and magnetic stress terms 
dominate the 5 term at late times and produce a slower 
decay of the shear distortion. Note that the presence of 
curvature anisotropy or any anisotropic trace-free mat¬ 
ter stress changes the shear evolution from the simple 
delta term that is usually studied in the literature (eg in 
Q). This behaviour can be seen in the dust era of a 
magnetic universe in the exact magnetic dust solution of 
Thorne and also in the study of almost isotropic dust 
plus collisionless trace-free radiation of HI- It was em¬ 
ployed in the study of nucleosynthesis in ref. in con¬ 
junction with the early microwave background isotropy 
data. Also, note that the ratio of the magnetic and black¬ 
body radiation densities is not constant (as assumed for 
example in and |16|), but falls logarithmically dur¬ 
ing the radiation era, and the three expansion scale- 
factors, a(t), 6 (t), c(t) evolve with time, to first order in 
the anisotropy, with a{t) oc h{t) oc t^/^{ln(t/to)}^ and 
c{t) oc t^/^{ln(t/to)}~^, with H = {2t)~^ as assumed 
at leading order. During the dust era we have a(t) cx 
h{t) oc — Kt~'^/^} and c{t) oc -|- 2Kt~'^/^}, 


with AT > 0 a constant such that ps/Pd = and 

H = 2/(3t) first order. 

The magnetic field and the accompanying shear dis¬ 
tort the microwave background temperature isotropy in 
accord with the most general anisotropic universes of 
Bianchi type VH. By combining eqn.(U) and the expres¬ 
sion for the angular anisotropy pattern in type VH, 0- 


ATA(f) = i^)isY (0, (/, Ho, / 150 , X, zis) (5) 
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where ‘Is’ denotes the epoch of last scattering {1 + zis = 
1100 with no reheating but if there is reionization of 
the universe then last scattering can be delayed until 
1 + zis = 78(flh/i5o)“^ for U.QZIS « 1 , M) where flf, is 
the baryon density parameter. The exact form of the pat¬ 
tern function Y can be read off from Eqs. (4.11)-(4.16) 
of Barrow, Juszkiewicz and Sonoda |l^]. The constant 
parameter x, introduced by Collins and Hawking |p4| is 
a measure of the three-curvature anisotropy configura¬ 
tion in the type VII universes. It contributes a spiral 
with 2/ttx twists to the angular pattern, ||^. If Hq < 1 
there will be a focusing the quadrupole towards the axis 
of anisotropy, generating a hot spot. 

The problem of constraining global anisotropy is sub¬ 
stantially different from the traditional statistical task of 
estimating parameters in Gaussian models. In the lat¬ 
ter case, the ensemble is entirely characterized by the 
power spectrum while in the former, a given set of pa¬ 
rameters corresponds to a completely specified pattern 
in the sky, up to an arbitrary rotation. This problem 
was dealt with in some detail in [Q. One can model 
the microwave background signal as the sum of two com¬ 
ponents: a statistically isotropic Gaussian random field 
AT/, which we assumed to have a scale invariant power 
spectrum on the scales we are interested in, and a global, 
anisotropic pattern, AT^, as in eq. (||), which is uniquely 
defined by the set of parameters x, flo, {a/H)o, and 9, 
4> (its orientation on the sky). Each pixel of the data set 
is given by di = (AT * (3){ri) + Ni where f3 represents 
the DMR beam pattern, is a unit vector pointing in 
the direction of pixel i and Ni is the noise in pixel i. To 
an excellent approximation, one can assume that N is 
Gaussian “white” noise, i.e. (NiNj) = sl5ij. 

Our task is, given a pair {x, Hq), to find the orien¬ 
tation ( 0 , (/>) which allows the maximum observed value 
of {a/H)o. One can do this using standard frequentist 
statistical methods: we define a goodness-of-fit statistic 
that depends on the data, compute its value for the actual 
data, and then compute the probability that a random 
data set would have given a value as good as the actual 
data. In [Q ry was defined to be: 


7] = min 771 ; 

<7,9,(p 


where rji = 


Al-Al 


( 6 ) 


Ag is the noise-weighted mean-square value of the data 
and Af is the noise-weighted mean-square value of the 
residuals after we have subtracted off the anisotropic 
part. Note that removing the incorrect anisotropic por¬ 
tion will only increase the residuals so the difference be¬ 
tween the two terms is an obvious choice for a goodness- 
of-fit. Dividing by Aq ensures a weak dependence on the 
amplitude of the isotropic component, while defining the 
statistic as the minimum of 771 allows us to deal with the 
uncertainty in ( 0 , </>). 

This statistical method was applied to the 4-yr GOBE 
DMR data-set. The two 53 GHz and the two 90 GHz 


maps were averaged together, each pixel weighted by the 
inverse square of the noise level, to reduce the noise level 
in the average map. All pixels within the Galactic cut 
were removed so as to reduce Galactic contamination, 
and a best-fit monopole and dipole were subtracted out. 
The map was degraded from pixelization 6 to pixelization 
5 (i.e. binning pixels in groups of four). Simulations were 
performed for a set of models from the (Hoj x) plane; for 
each choice of the three parameters (Hq, x and {a/H)o) 
approximately 200 to 500 random DMR sets were gener¬ 
ated, so allowing us to determine an approximate fit to 
the probability distribution function of 77 . 

The most conservative limit on the cosmological mag¬ 
netic field arises when we assume pga < pb so the whole 
anisotropy is contributed by the magnetic field stresses. 
A reasonable fit of the upper bound at a 95% confidence 
level is: 

( —)o<^^^/(x,Ho)x10-® 
p 4 

/ ~ 2.1 

with X e [.01, 3] and Hq S [.1,1], (7) 

where we have considered the largest possible contri¬ 
bution from the magnetic component. Note that the 
“shape” factor is roughly bounded by 0.6 </ < 2.2. This 
gives us a final bound on the magnitude of the magnetic 
field today of 

Bq < 3.5 X 10 “®/i(Ho/ 75 Q)iGauss. ( 8 ) 

This bound can be improved by a factor of -s/S if one 
considers the results from Q. In this case, a slightly 
different goodness-of-fit statistic is used: instead of work¬ 
ing with the noise-weighted quantities, Ag and Af, the 
authors chose to weight the pixels with the covariance 
matrix of the total Gaussian components (i.e. the noise 
and isotropic cosmological components). This gives a 
limit of 

Bq < 2.3 X 10“®/5 (Ho/75g)5Gauss (9) 

Note that the microwave background limits on the am¬ 
plitude of anisotropies are much stronger than those im¬ 
posed by nucleosynthesis Q . In unrealistic models with 
no anisotropic matter stresses (and therefore no magnetic 
field) and isotropic curvature, the shear falls off rapidly 
in accord with the d term in (^) and the limits from 
nucleosynthesis would be stronger. But with anisotropic 
matter stresses, magnetic fields, or anisotropic curvature, 
the anisotropy falls only logarithmically during the radia¬ 
tion era. The limits on pb/P r at nucleosynthesis are only 
0 ( 1 ) and the log decay means they are weaker than limits 
0 ( 10 “®) imposed at z = 1.1 x 10 ® by the microwave back¬ 
ground. If there is reheating and last scattering occurs at 
2 ; << 1100 then the analysis is slightly changed but last 
scattering would need to occur at a redshift lower than 
6 Ho ^50 for ffio nucleosynthesis limit to be competitive 
with the microwave limit. This never happens. 
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Adams et al have used the weak nucleosynthesis 
limit of Grasso and Rubinstein of Bq < 3 x 10“^ Gauss 
obtained for random fields to argue that a cosmologi¬ 
cal magnetic field could lead to observable distortions of 
the acoustic peaks in the microwave background. Our 
limit on Bq rules out any observable effect on the acous¬ 
tic peaks. In fact, when the (lnt)“^ decay of pbIpt of 
eqn.(^ is taken into account a nucleosynthesis limit is 
strong enough to render the acoustic peak distortions 
unobservable. Our limit permits a field strength of 10“® 
Gauss required to induce a measurable Faraday rotation 
in the polarization of the microwave background ||^ . 

Any period of inflation long enough to explain the 
horizon and flatness problems would necessarily reduce 
magnetic field effects and their associated anisotropies 
to unobservably low levels. If N e-folds of de Sitter in¬ 
flation occurs {p = —p =constant) then a/H will be 
reduced by exp(—3A^) and ps/p will be reduced by a 
factor exp(—4A^) and A^ ^ 70 is sufficient to solve the 
horizon problem |^. The formula for a/H ap¬ 
plies to the case of inflation if the 5 term is changed 
to 5exp(—3A^). Note that if generalized inflation occurs 
(0 < 7 < 2/3), the 5 term of eqn. (^) always domi¬ 
nates the {pga pb)/P term ast^oo unlike in the non¬ 
inflationary case when 2/3 < 7 < 4/3. All anisotropies 
decay in accord with the no-hair theorem when the cur¬ 
vature is non-positive because the anisotropic trace-free 
stresses obey the strong energy condition |^. The dis¬ 
covery of microwave background patterns characteristic 
of large-scale homogeneous anisotropy or a homogeneous 
primordial magnetic fields in future observational pro¬ 
grams would therefore rule out the standard picture of 
inflation. 
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